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O ABSTRACT 



We present a far-infrared galaxy luminosity function for the local universe. We have 



obtained 389 spectroscopic redshifts for galaxies observed at 90 /xm in the AKARI 
i ^ i Deep Field South, using the AAOmega fibre spectrograph via optical identifications 

in the digitized sky survey and 4m-class optical imaging. For the luminosity function 
i-H presented in this paper, we have used those galaxies which have redshifts < z < 0.25, 

£> have optical magnitudes and are not part of a newly discovered cluster of galaxies 

<3\ (giving a total of 130 sources). Infrared and optical completeness functions were es- 

timated using earlier Spitzer data and APM B-band optical data respectively, and 
the luminosity function has been prepared using the 1/Vmax method. We also sepa- 
rate the luminosity function between galaxies which show evidence of predominantly 
star-forming activity and predominantly active galactic nucleus (AGN) activity in 
their optical spectra. Our luminosity function is in good agreement with the previ- 
ous 90 /zm luminosity function from the European Large Area ISO Survey, and we 
also present a luminosity function with combined AKARI and ISO data. The result 
is in reasonable agreement with predictions based on the earlier IRAS 60 ^m PSC-z 
catalogue, and with a recent backward evolution model. 
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1 INTRODUCTION 



Recent advances in observational infrared astronomy have 
led to a rapid increase in the constraints which can now be 
placed on theoretical models of the development of struc- 
ture in the Universe. The discovery of far-infrared luminous 
populations of galaxies by SCUBA (e.g. Hughes et al. 1998, 
Smail et al. 1997) and of dust obscured galaxies by Spitzer 
(Dey et al. 2008) have led to the realisation that massive 
galaxies formed the bulk of their stars at high redshifts dur- 
ing bursts of star-formation, contrary to the simpler growth 
expectations predicted by earlier hierarchical models (see 
review in Hauser & Dwek 2001). Furthermore, there has 
been a growing realisation that feedback from intermittent 



AGN activity is important in regulating stellar mass assem- 
bly (e.g. di Matteo et al. 2005). 

The evolution of star-forming galaxies (SFGs) and 
galaxies with active galactic nuclei (AGNs) can provide im- 
portant observational constraints on galaxy evolution mod- 
els. Radiation from AGN dust tori is predicted to peak in the 
mid-infrared (Efstathiou & Rowan- Robinson 1995), whereas 
dust radiation from the giant molecular clouds which domi- 
nate the bolometric output from star formation means that 
SFGs will peak in the far-infrared, between 50 fim - 200 /im. 

There has been considerable work recently in the in- 
frared shortward of this SFG peak. A number of key pa- 
pers on the history of star formation and the evolution of 
the luminosity function have used 24 /im data from the 
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Spitzer Space Telescope. An extensive study of the evolu- 
tion of star-forming galaxies showed an order of magnitude 
increase in typical infrared luminosity from the present out 
to z — 2 (Perez-Gonzalez et al. 2005). A recent study of 
Spitzer 8 (im - 24 /im sources suggested a rapid increase of 
the co-moving volume emissivity back to z ~ 1 and a con- 
stant average emissivity at z > 1 as well as showing evidence 
of downsizing in star formation (Rodighiero et al. 2010). 
Magnelli et al. (2009) presented luminosity functions at 15 
/im (based on Spitzer 24 /im data) and at 35 /im (based on 
Spitzer 70 /im data) between redshifts 0.4 < z < 1.3. Ru- 
jopakarn et al. (2010) presented a 24 /im luminosity func- 
tion in the range < z < 0.65, and also combined their 
results with those of Magnelli et al. (2009) to cover red- 
shifts out to z — 1.2, confirming strong evolution, consis- 
tent with pure luminosity evolution. Several studies have 
found that star formation from Luminous InfraRed Galax- 
ies (LIRGs, lO 11 L <L / j i <lO 12 L0) and Ultra-Luminous 
InfraRed Galaxies (ULIRGs, 10 12 L Q <L /ii <10 13 L Q ) have 
made a much larger fractional contribution to the co-moving 
volume-averaged star formation density at z = 1 than they 
do locally (e.g. Le Floc'h et al. 2005, Caputi et al. 2007). 
A study of data from the AKARI Space Telescope's Deep 
Field North based on observed wavelengths of 2 /im - 24 /im 
which analysed data out to z ~ 2.2 has also found strong 
evidence for evolution of galaxy luminosities towards higher 
redshift (Goto et al. 2010). 

Work has also been done recently at wavelengths longer 
than the SFG peak. A measurement of the ISO 170 /im lumi- 
nosity function was made by Takeuchi et al. 2006. Luminos- 
ity functions based on observations at 250, 350 and 500 /mi 
from the recent Balloon-borne Large Aperture Submillime- 
ter Telescope (BLAST) survey centred on the GOODS- 
South field have shown a strong increase in the number 
density of the most luminous galaxies from the present out 
to z = 1 (Eales et al. 2009). Two recent papers have pre- 
sented luminosity functions using observations at 250 /im 
from the science demonstration phase of the Herschel Space 
Telescope: Dye et al. (2010) used data from a 14 deg 2 field 
out to z = 0.5, showing steady evolution out to this redshift, 
and Eales et al. (2010) in two smaller but deeper fields to 
z — 2.0 found strong evolution out to z ~ 1 but at most 
weak evolution beyond this redshift. 

The first measurements at wavelengths closer to the 
SED peaks of local star-forming galaxies were made by 
the Infra-Red Astronomical Satellite (IRAS). Saunders et 
al. (1990) presented a local luminosity function of 60 fim- 
selected IRAS sources and found evidence of strong evolu- 
tion (which could be explained by evolution in either lu- 
minosity or density). This work was extended and revised 
by Soifer & Neugebauer (1991), Sanders et al. (2003) and 
Takeuchi et al. (2003). Wang & Rowan- Robinson (2010) 
confirmed this evolution using a much larger sample of 
IRAS 60 /jm sources in the redshift range 0.02 < z < 0.1. 
The evidence from the ISO ELAIS 90 (im survey suggested 
a calibration disagreement with the IRAS 100 /im data 
(Heraudeau et al. 2004) so it is important to obtain in- 
dependent estimates of the local luminosity function, par- 
ticularly at 90 - 100 /im. A 90 /im luminosity function was 
derived from the European Large Area ISO Survey (ELAIS 
paper IX, Serjeant et al. 2004), based on spectroscopic red- 
shifts of 151 Infrared Space Observatory (ISO) sources in 



the Northern ELAIS fields. With optical identification of 29 
ISO far-infrared (FIR) sources in the Lockman Hole using 
the FIR-radio correlation for positional accuracy, Oyabu et 
al. (2005) showed that the space density of the FIR sources 
at 0.3 < z < 0.6 was some 460 times larger than the lo- 
cal value. A recent luminosity function using Spitzer 70 jim 
sources in the ultra-deep GOODS-North pencil-beam survey 
showed evidence for evolution at redshifts z > 0.4 (Huynh 
et al. 2007). An analysis to find the split between SFGs 
and AGN-powered galaxies in Spitzer 70 /im data in the Ex- 
tended Groth Strip field found that virtually all were pow- 
ered predominantly by star-forming activity (Symeonidis et 
al. 2010). 

AKARI's Far Infrared Surveyor Instrument (FIS) is 
most sensitive at 90 /jm, and provides the opportunity to ex- 
pand the observational evidence in this region near the SFG 
peak. This paper aims to provide a new 90 /im luminosity 
function, to compare with the previous ISO one, and also 
to combine with data from that survey to provide the best 
luminosity function at this wavelength so far. The present 
study aims to measure a local benchmark against which to 
compare luminosity functions of more distant SFGs based 
on more detailed, deeper data which will become available 
in the near future from instruments such as Herschel Space 
Telescope. In this paper, we will constrain the 90 /im galaxy 
luminosity function for z < 0.25, using the spectroscopic 
redshifts we have found with the AAOmega fibre spectro- 
scope for sources in the AKARI Deep Field South (ADF-S, 
Matsuura et al., 2011). An earlier version of this luminos- 
ity function (using about one-half of the data now available, 
showing similar results) was included as part of a confer- 
ence proceedings (Sedgwick et al. 2009). In a separate paper 
(Goto et al. 2011) we are estimating the local far infrared 
luminosity function by cross-correlating the AKARI all-sky 
survey with IRAS Bright Source catalogue. 

The methods used to find the infrared sources, to iden- 
tify redshifts from cross-identified optical spectra of these 
sources and to derive the luminosity function, are described 
in Section [2] The results are presented in Section [3] and dis- 
cussed in Section ID 

This paper assumes the parameter values of the current 
concordance cosmology: _ffo=72.0 km s~ Mpc^ 1 , Qm~0.3 
and Qa=0-7- 



2 METHODS 

2.1 AKARI Deep Field South Observations 

The AKARI Space Telescope (Murakami et al. 2007) is a 
cryogenically cooled 68.5 cm mid/far-infrared telescope in 
a low Earth, Sun-synchronous orbit. Besides performing an 
all-sky survey at 10 /xm-160 /im, AKARI also conducted a 
series of deep pointed observations. These deep guaranteed 
time extragalactic survey programmes were conducted near 
the North and South Ecliptic Poles for reasons of orbital 
visibility, and are described in Matsuhara et al. (2006) . Both 
the FIS (60 - 160 /im) and the InfraRed Camera (IRC, 2 - 
24 /im) instruments were used in these surveys. 

The AKARI Deep Field South (ADF-S) is centered on 
RA 4h 44m 00s, Dec -53° 20' 00" (J2000), a low-cirrus 
region near the South Ecliptic Pole. The FIS survey in this 
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Table 1. Cross-identifications between the AKARI FIS 90 fim 
catalogue and other infrared catalogues in the ADF-S. The la 
positional uncertainty shown below is based on the formula 
5RA=<5DEC=0.6(1/SNR)FWHM (Ivison et al. 2007) for FIS and 
Spitzer. The BLAST uncertainty is quoted in Valiante et al. 
(2010). Therefore a 2<r radius for a combination of sources is be- 
tween about 10" and 15", so matches at these distances are shown 
below. For AKARI FIS details (data not yet published), see Shi- 
rahata et al. (2009). For BLAST data and results, see Valiante et 
al. (2010). For Spitzer results (data not yet published), see Scott 
et al. (2010). AKARI IRC data has not yet been published. 



Table 2. AAOmega observations: exposure times and number of 
redshifts identified from emission lines. 



Date Total Exposure Redshifts 





hr:m:s 


Identified 


29 Oct 2007 


2:30:00 


130 


30 Nov 2007 


1:23:40 


63 


3 Jan 2008 


2:00:00 


59 


27 Nov 08 a 


1:46:40 


76 


27 Nov 08 b 


0:26:40 


61 


Total 




389 



Matches to 
FIS 90 /im 





Total 


lcr 


sources 


within 


Catalogue 


Sources 


Error 


10" 


15" 


AKARI FIS 90 /jm 


2,282 


5.1" 






AKARI FIS 65 fim 


391 


5.4" 


125 


171 


AKARI FIS 140 fim 


315 


7.0" 


56 


80 


AKARI FIS 160 fim 


216 


6.5" 


7 


11 


Spitzer 70 fim 


1,205 


5.2" 


583 


72 


Spitzer 24 fim 


41,503 


1.1" 


1,199 


1,584 


AKARI IRC 15 fim 


725 


2.0" 


25 


28 


AKARI IRC 24 ^m 


184 


2.0" 


24 


30 


BLAST 250/350/500 /im 


232 


5.0" 


22 


38 



region covered over 12deg 2 .IRC imaging was also conducted 
in Open Time over part of this field (about 0.8 deg 2 ). 

The sources reported in this paper were detected in the 
ADF-S by the 90 /im WIDE-S band of the FIS, using its 
slow-scan observation mode with point-source photometry 
(Shirahata et al. 2009). This reached 28mJy at 60% com- 
pleteness (our infrared completeness function is based on 
earlier Spitzer data: see Section 2.6 for details). By way of 
comparison, a 50% completeness limit of the sources used 
in the 90 fim luminosity function for the ELAIS Northern 
fields reached 70 mjy (Serjeant et al. 2004). 

The FIS survey found a total of 2,282 sources at 90 fim 
at a signal-to-noise ratio SNR>5 (giving a minimum flux 
of 12.81 mjy) in this field (Shirahata et al., in prepara- 
tion). Cross-identifications with other infrared catalogues 
are shown in Table [l] 



2.2 Optical Identifications 

We used the Bj-band APM catalogue for optical cross- 
matching of FIS sources, and for estimating the optical com- 



pleteness of the FIS catalogue (discussed in Section 2.6 1. 
This allows a comparison with earlier ELAIS work, which is 
in the Bj-band. 

Optical identifications of the AKARI sources were 
checked by eyeballing all the candidate sources using R- 
band images from the Digitised Sky Survey (DSS) and 
(where available) images from the CTIO 4m telescope. We 
chose DSS-R for eyeballing identifications of both IRC (mid- 
infrared) and FIS sources, partly because this better suited 
the IRC population and partly to match the available CTIO 
R-band data over part of the field. Where available, identifi- 
cations were supported by sub-arcsecond positions from the 
1.4 GHz Australia Telescope Compact Array survey (White 



et al. in preparation) and the AKARI-IRC data (Pearson 
et al. in preparation). Since we are dealing in this paper 
with bright local FIS (far-infrared) galaxies, identifications 
in DSS B-band have counterparts in the DSS-R band, and 
vice versa. 



2.3 AAOmega Spectroscopic Redshift 
Identifications 

We used AAOmega, the fibre-fed optical spectrograph at 
the Anglo Australian Observatory, to obtain optical spectra 
of a selection of the ADF-S sources in order to determine 
spectroscopic redshifts. For our targets, we used a sparse 
sample based on the AAOmega fibre configuration, which 
(in the absence of any evidence to the contrary) we have 
assumed to be unbiased. The area covered was n degrees 
(one AAOmega field-of-view: see Figure Only sources 
for which we had optical counterparts were targeted, and 
these were cross-identified with the AKARI source coordi- 
nates by eye. A few very local sources were extended op- 
tical objects, and the cross-identified AKARI coordinates 
were within the extended optical image. Most of the cross- 
identified optical sources were point sources, 94% of which 
fell within 10", and all fell within 16"of the AKARI coor- 
dinates. The AKARI telescope pointing error is less than 
3" (Murakami et al. 2007). 

Data were obtained by AAOmega during four sessions 
between October 2007 and November 2008. None of the 
nights was photometric. Gratings were 385R and 580V with 
central wavelengths 7250 A and 4800 A, and wavelength res- 
olution A/AA=1300. The spectra were reduced using the 
standard AAOmega pipeline, which was also used to com- 
bine the red and blue arms. The combined spectra each have 
a range of about 3800 - 8800A, although noise at either end 
of this range, particularly the blue end, often reduced the 
useful data to a range of about 4000 - 8700A. 

Redshifts were measured using the following emission 
lines: [Nil] /Ha/ [Nil] AA6548, 6563, 6583 and the [SB] dou- 
blet AA6716, 6731, and the three H/3/[OIII]/[OIII] AA4861, 
4958, 5007 emission lines. Redshifts were only taken where 
one or other (or both) of these two triplets of emission lines 
were identified. Unless it was redshifted out of range, the 
[SII] pair was invariably seen with the [Nil] /Ha/ [Nil] triplet. 
In addition, the asymmetrical [OH] emission line resulting 
from the close doublet AA 3726, 3729 was often found at 
the same redshift to support the identifications. Examples 
of spectra obtained showing the redshift identifications are 
given in Figure [l] 
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AAOmsda Nov-OB run 1, Fibre Number 324 




FIS Redshift Sources in the AKARI Deep Field South 



5000 6000 
Wavelength / Angstroms 

AAOmega Nov— OB run 1, Fibre Number 155 



6000 7000 
Wavelength / Angstroms 



Figure 1. Examples of spectra which show the main emis- 
sion lines used to identify redshifts in this paper: from left, 
[OII]A3726, H/3/[OIII]/[OIII] AA4861, 4958, 5007, [Nil] /Ha/ [Nil] 
AA6548, 6563, 6583 and [SII]AA6716, 6731 (rest-frame wave- 
lengths). Note that the strong emission line at ~5570A is an 
artefact from street lighting in a local town; there are also some 
telluric lines (e.g. at ~7700A). The spectra have been normalised 
to 1 using boxcar median smoothing to remove the effect of dif- 
fering throughputs by the blue and red arms of the spectrograph. 
Using the Kauffmann criterion on the BPT diagram as described 
in the text, the top source (which is at z = 0.071) is an SFG, and 
the bottom source (z = 0.211) is an AGN. 



The quality of the spectra obtained varied widely be- 
tween the different nights, as a result of varying asmospheric 
conditions and length of exposure times achieved for each set 
of targets (see Table[2j. Overall, a total of 389 redshifts were 
measured for the 1179 distinct sources targeted. Some of 
these were for non-FIS sources (e.g. AKARI-IRC, BLAST) 
and will be included in a future paper. For the FIS sources 
used in this study, a total of 317 sources for which we ob- 
tained optical counterparts were targeted, and redshifts were 
obtained for 183, a success rate of about 58%. Of these, 152 
had redshifts in the range < z < 0.25 and also had optical 
magnitudes, and these (subject to the removal of a new clus- 



ter described in 2.5 1 will be used in preparing the luminosity 
function presented in this paper, and are shown in Figure [2] 
Sources without spectroscopic redshifts were not used. 

The distribution of luminosity against spectroscopic 
redshifts used in preparing this luminosity function is shown 
in Figure [3] (which also highlights the new cluster found at 
z ~ 0.06 discussed in Section 2.5 below). 



2.4 AGN / SFG Discrimination 

Two indicators were used to discriminate between AGNs 
and SFGs in our sources. 

(a) Quasars are often defined as sources which have 




4.a 4.6 

FtA / hours 

Figure 2. FIS sources with z<0.25 (red plus symbols) on all FIS 
sources (black dots) in the AKARI Deep Field South. 



Luminosity v Redshift (with cluster in red) 




0.10 0.15 
Redshift 



Figure 3. The redshift - luminosity plane for all sources used, 
plus the cluster sources (which arc highlighted in red). Luminosity 
is i/L^/Lq as calculated for the luminosity functions in this paper. 



broad Balmer emission lines with FWHM over 1000 km s~ . 
Although we were able to measure FWHM of the Ha line 
for all the 152 selected sources, only three of them had such 
broad emission lines. In this paper, we follow Vanden Berk 
et al. (2001), and regard sources with FWHM of the Ha line 
over 500 km s" 1 as indicating significant AGN activity. A 
total of 10 of the sources in our sample showed this. 

(b) A BPT diagram (Baldwin, Phillips & Terlevich 
1981) is plotted in Figure[4] relating the two intensity ratios 
of [Nil] A6583 / Ha A6563 and of [OBI] A5007 / U/3 A4861. 
The black line discriminates between AGN sources and SFG 
sources, based on the formula in Kauffmann et al. (2003), 
with a position above this line indicating an AGN. A total 
of 105 sources had sufficient spectral data to plot the ratios 
in this diagram, of which 39 were identified as AGNs by the 
Kauffmann line. 

Most of the sources classified as AGNs were therefore 
selected by the BPT method, although many of these were 
not broad-line AGN since they had FWHM of the Ha line 
< 500 km s _1 . In addition, there were 5 sources with the 
Ha line > 500 km s , 3 of which did not have the ratios to 
use the BPT diagram, and 2 of which had ratios below the 
BPT Kauffmann classifying line. Classifying these additional 
sources as AGNs did not significantly affect the separate 
luminosity functions in Section [3] below. 

The ratio of [OIII] A5007 /[Oil] A3727 was also consid- 
ered as a possible AGN/SFG discriminator. However, this 
ratio proved unreliable: the two emission lines are not close 
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BPT Diagram of selected AKARI-DFS sources 
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Figure 4. BPT diagram of the 105 sources used in the luminosity 
functions which had data to calculate these two ratios. The black 
line plots the Kauffmann et al. (2003) formula to discriminate 
between AGNs and SFGs. The area between this black line and 
the red line contains AGN-SFG composite sources, based on the 
formulae in Kewley et al. (2006). 
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Figure 5. The new cluster found in the survey field. Cluster 
members are shown as filled red diamonds; sources in the same 
redshift range but not cluster members as filled blue diamonds; 
other redshift sources as open diamonds. 



and are from the two (blue and red) arms of the spectro- 
graph which have different throughputs; in addition, there 
could be reddening intrinsic to the galaxies. In fact, this in- 
dicator usually conflicted sharply with the other indicators. 
In contrast, the ratios used in the BPT diagram are of lines 
which are close in wavelength, so are not affected either by 
throughput variations or by reddening. 

The result of this analysis gave 44 optically-selected 
AGN sources and 108 optically-selected SFG sources out of 
the total of 152 FIS sources with z < 0.25 and optical mag- 
nitudes. Note, however, that many of our sources lie close 
to the dividing line, and the Kewley et al. (2006) criterion 
suggests many are AGN-SFG composite sources (Figure Ek. 
Both star formation and active nuclei contribute to the bolo- 
metric energy output of these galaxies, although the BPT 
diagram on its own is insufficient to determine bolometric 
fractions. 



2.5 Cluster of Galaxies in the Field 

As can be seen in the luminosity / redshift plot (Figure |3k, 
there is a group of sources bunched at about z = 0.06. In 
fact, there are 22 sources at a redshift of between 0.056 < 
z < 0.061 and a separation on the sky of under 45 ar- 



cminutes, centred on (04 42 12.5, -53 30 46). This redshift 
range corresponds to a depth of 21 Mpc, and this angular 
separation corresponds to a width of about 3.2 Mpc. The 
cluster sources are identified in Figures [3] and [5] and listed 
m Table [3] Figure [3] shows that the cluster contains rela- 
tively more luminous sources than the sample as a whole. 
This may have been caused by under-representation of high- 
luminosity galaxies at the lowest redshifts in our sample, or 
we may have sampled a region of enhanced star-formation 
in the outer (infalling) region of the cluster (Coppin et al. 
2011, but see Bai et al. 2009 and Haines et al. 2011). Includ- 
ing this cluster in the sample distorts the luminosity function 
and causes the mean of V/Vmax to fall well below 0.5 (this is 



discussed in Section 3.2 below). Therefore, we have excluded 
the sources in the cluster from the final luminosity function. 
(This is not the well-known Dressier DC0428-53 cluster at 
z — 0.041, which is within the ADF-S field but outside the 
part of the field for which we have obtained redshifts.) 

We have estimated the comoving volume of the cluster 
as ~10 4 Mpc 3 . For this calculation, we converted the red- 
shift range to comoving depth per unit area and multiplied 
by the area of the cluster. This cluster volume is a tiny 
proportion (~ 0.01) of the whole field, so excluding the 
cluster volume from the calculations makes no significant 
difference to either the luminosity function or (V/Vmax) so 
we have not made this adjustment. 



2.6 Luminosity Function Methods 

The 90 (jm local luminosity function was estimated as fol- 
lows. 

Completeness as a function of flux based directly on 
the FIS data is still under construction (Shirahata et al., in 
preparation). Therefore we have had to make our own as- 
sessment of the completeness of the AKARI data. The ultra- 
deep Spitzer 70 (jm observations of the GOODS-North field 
reached down to a flux of 1.2 mjy (Frayer et al. 2006). This 
is an order of magnitude lower than the limit of the present 
AKARI survey, so we assume this will be 100% complete 
down to the lower end of the AKARI data. Frayer et al. 
(2006) showed that this Spitzer data is well fitted by the 
model of Lagache et al. (2004) which we used to generate 
number counts against which to measure the AKARI com- 
pleteness. However, since this data is at 70 fim we have had 
to make an adjustment for our 90 (im case, and to do this 
we used the model described in Pearson (2001), scaling the 
fluxes by a factor of 0.673 on the basis of that model's pre- 
diction that 7V(SVo M m > 0.0673 Jy) = N(S 90 ^ m > OTJy). 
The resulting completeness ranges from 97% at 100 mjy to 
60% at 28 mjy, as shown in Figure [6] 

As mentioned earlier, sources without APM B-band 
counterparts were excluded. Using APM B-magnitudes, we 
estimated spectroscopic completeness as a function of mag- 
nitude, using the AAOmega redshift sources compared to 
APM sources cross-identified with the sources in the total 
ADF-S FIS sample before identification of redshifts. Each 
point in Figure [7] shows the completeness using a range of 
±1.5 magnitudes. The lcr error bars of the points are shown 
in the figure. The completeness is close to 100% down to 
B-band magnitude of about 18 after which it declines to 0.1 
at about magnitude 21. 
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Figure 6. Estimated FIS completeness function using Lagache 
model (Lagache et al. 2004) fitted to Spitzer 70 fim GOODS-N 
data, modified for the 90 (im case by the Pearson model (Pearson 
2001) as described in the text. 
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Figure 7. Estimated AAOmega optical completeness function 
using APM catalogue. Vertical lines show la errors, and horizon- 
tal lines show bin sizes. 



For galaxies with B>21 and z<0.25, we have adopted 
the assumption used for ELAIS (see appendix to Serjeant 
et al. 2001 for a detailed discussion) that such galaxies have 
more local counterparts, and at least one member of such a 
class of galaxy would be included in our sample. This would 
mean that its contribution to the luminosity function using 
the 1 /Vmax method would be given a weighting that exactly 
accounts for the increased difficulty of detecting such a class 
of galaxy. 

To check whether or not our sample is typical of the 
sky as a whole - i.e. whether it is affected by large scale 
cosmic variance - we compared the number per unit area 
with that obtained in the 50 deg 2 Spitzer Wide-area Infrared 
Extragalactic Legacy Survey (SWIRE) at 70 (jm (quoted 
in Frayer et al. 2009) for > O.Uy objects (following the 
procedure adopted in ELAIS, e.g. Serjeant et al. 2001). The 
fraction of number predicted / actual number was 1.80±0.32 
and we renormalised the effective luminosity density (and 
error) by this fraction. 

As explained above, we have ended up with a final total 
of 130 sources (after excluding the cluster) for the luminosity 
function. For comparison, the fields used in the ELAIS IX 
luminosity function had a total of 151 sources (as with our 
sample, these all had spectroscopic redshifts). More sources 
would obviously have been welcome, and we hope to obtain 
further AAOmega observations in the future. This survey 
covers a similar redshift range to the ELAIS IX survey (cf. 
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Figure 8. Local galaxy luminosity function using all ADF-S 
sources with AAOmega spectroscopic redshifts selected as de- 
scribed in the text. It is assumed in this figure that there is no 
pure luminosity evolution. The hatched area shows Poisson error 
bands. The line is the double power-law prediction by Serjeant & 
Harrison 2005. 



Figure 1 of Serjeant et al. 2004), so in addition to showing 
our new results (Figures [8] to 10 1, we have felt justified in 
combining the two samples. The combined luminosity func- 
tion is shown in Figure [IT] and discussed below. 

We calculated the 1/Vmax luminosity function (Schmidt 
1968) following the methodology in Serjeant et al. (2004). 
K-corrections were made assuming the M82 star-forming 
spectral energy distribution, except for the AGN luminosity 
function for which we assumed a Seyfert 2 SED template 
derived from the Spitzer Wide-Area Infrared Extragalactic 
Survey (SWIRE). 



3 RESULTS 

3.1 ADF-S Results 

As expected from earlier studies (Le Floc'h et al. 2005, Ca- 
puti et al. 2007), our z < 0.25 sample is largely free of LIRGs 
and ULIRGs. The highest luminosity was 1.3 x 10 11 L , and 
only two sources were over 10 11 Lq, the usual low-end def- 
inition of LIRGs (Sanders & Mirabel 1996), where we use 
idj„(90^im) as a rough estimate of the bolometric infrared 
luminosity. 

The local 90 fim luminosity function for all sources is 
shown in Figure [8] and tabulated in Table [4] For each lumi- 
nosity bin, we have calculated the space density of objects 
with luminosities within ±0.15 dex. 

The red line in Figures [8] to [To] is the double power- 
law prediction (Equation [T| for the local 90 /im luminosity 
function from Serjeant & Harrison (2005), which was based 
on extrapolated fluxes of 15,411 local galaxies (z < 0.1) 
in the PSCz catalogue (Saunders et al. 2000) derived using 
IRAS data: 
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Figure 9. Star-forming galaxy luminosity function (with AGNs 
excluded). The line is the double power-law prediction by Serjeant 
& Harrison 2005 (for all galaxies). 



3.2 Mean of V/Vmax 

The mean of V/Vmax should be equal to 0.5 if the sample is 
complete, assuming a uniform space density. We found that 
(V/Vmax) = 0.478±0.022 if we assume that there is no pure 
luminosity evolution. However, if (1+z) 3 pure luminosity 
evolution is assumed, (V/Vmax) falls sharply to 0.427±0.023, 
although the luminosity function itself is not markedly dif- 
ferent (compare Figures [8] and 11 left). Much of this short- 
fall is accounted for by the sources at 0.2< z <0.25: making 
a cut at z=0.2 increases (V/Vmax) from 0.427 to 0.467 for 
the luminosity function with pure luminosity evolution of 
(1+z) 3 . (This reduces the number of sources of sources from 
130 to 119.) 

We have used this value for pure luminosity evolution 
because earlier work has usually found pure luminosity evo- 
lution of this order, albeit with large uncertainties: for ex- 
ample, ELAIS IV found (l+z) 2 ' 45±a85 (Serjeant et al. 2001) 
and ELAIS IX found (l+z) 3 ' 4±1 '° (Serjeant et al. 2004), 
both at 90 fxm, and Saunders et al. (1990) found (l+z) 3 0±10 
at 60 /jxa. 
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Figure 10. AGN luminosity function (with SFGs excluded) as 
described in the text. The line is the double power-law prediction 
by Serjeant & Harrison 2005 (for all galaxies). Note that only 36 
AGNs were identified in the sample of 130 sources. 



3.3 Combined ADF-S and ELAIS Results 

The 90 /im luminosity function for the Final Analysis of the 
northern ELAIS fields was presented in Serjeant et al. (2004) 
and is reproduced in the central panel of Figure [IT] This 
covered three Northern ELAIS fields (totalling 7.4 square 
degrees) to a flux limit of 70 mjy. Completeness reached 
100% at a flux of 150 mjy and 50% at 70 mjy, so it was 
not as deep as the ADF-S survey. The ELAIS survey was 
specifically designed to minimise the effects of cosmic vari- 
ance (Oliver et al. 2000) and so was spread over three fields, 
giving confidence that no cosmic variance renormalisation 
was necessary. 

The left panel in Figure [TT] shows the ADF-S luminosity 
function recalculated assuming pure luminosity evolution of 
(1 + z) 3 in order to compare and combine like-with-like with 
ELAIS, although as will be evident this makes only a slight 
difference to our luminosity function at this range (z < 0.25) 
of redshifts. The right-hand panel of this Figure shows the 
combined (ADF-S and ELAIS) luminosity function. 



where log 10 = — 2.93dex 1 Mpc 3 , the character- 
istic break luminosity L, is given by log 10 L» = 
23.54 WHz" 1 sr"\ /3 = 0.567 and 7 = 2.84. This prediction 
assumes pure luminosity evolution of (1 + z) 3 . Converting 
to the units used in this paper (and using the more current 
value H = 72) gives the predicted L, = 3.06 x 10 10 L . 

The luminosity functions for star-forming galaxies (af- 
ter excluding all those with any evidence of AGN activity) 
and for AGN sources are shown separately in Figures [9] and 
|10| In the case of AGN sources, the boxcar width is ±0.2 
dex due to the lower number of data points. 



4 DISCUSSION 

The most noticeable departure from the predicted local 90 
/xm luminosity function of Serjeant & Harrison (2005) de- 
scribed above is a deficit of galaxies from just below L* 
(from ~ 1Q w Lq). This is particularly clear in the luminosity 
function where galaxies with AGN have been excluded (Fig- 
ure |9J, where it is significant at about 3a. This may suggest 
we are sampling atypical large-scale structures, and provides 
strong motivation for expanding our spectroscopic area to 
overcome these effects. To put this in context, the comoving 
distance across the AAOmega field at a median redshift of 
z ~ 0.1 is about 14 Mpc. The ELAIS results (Figure [TTja) 
did not show the deficit of L„ galaxies which we have found 
in this study, and the combined luminosity function (Fig- 
ure |11| :) is close to the Serjeant & Harrison (2005) double 
power-law prediction. 

At the faint end of the luminosity function, there is a 
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Figure 11. Comparison between 90 fim luminosity functions for (from left) ADF-S, ELAIS and combined ADF-S and ELAIS. 
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Figure 12. Comparison of the AKARI ADF-S 90 fim luminosity function with predictions from the new Gruppioni et al. (2011) model. 
The AKARI data points are independent in this figure (unlike the boxcar-smoothed data in Figures |8)ll| . The plots from the Gruppioni 
et al. (2011) model are: black line = total, cyan= spiral galaxies, red= starburst galaxies, bluc=Scyfert 2s, green=obscured AGNs and 
brown = AGN1 sources. The left figure shows all AKARI sources, the centre figure shows SFGs compared with the two star-forming 
Gruppioni classes, and the figure on the right shows AGNs. 



slight excess number density when all sources are consid- 
ered compared to the Serjeant & Harrison (2005) prediction 
(see Figure [8}. At the bright end, we find that the slope is 
shallower for AGNs than for SFGs (cf. Figures [9] and [10} in- 
dicating an increased contribution at high luminosity from 
AGN sources, which we might expect given the inferences 
from mid-infrared spectroscopy and radiative transfer mod- 
els of increased AGN bolometric fractions with increasing 
bolometric luminosities (see for example Genzel et al. 1998 
and Genzel & Cesarsky 2000). The relative number density 
of AGN sources rises from under 10% of the total number 
density at the faint end to over 30% of the number density 
of all sources at the bright end. 

We have also compared our observations with the the- 
oretical prediction for the 90 /im luminosity function for 
< 2 < 0.25 galaxies from a new backward evolution model 
of Gruppioni et al. (2011) which uses the methodology devel- 
oped in Gruppioni et al. (2008). This comparison is shown 
in Figure [12] for all sources and for SFGs and AGNs sep- 
arately. The Gruppioni model includes separate evolution 
for five types of galaxy: spirals (with low to moderate star 
formation), starbursts, low-luminosity AGNs (Seyfert 2s), 
AGNls, and AGN2s (obscured AGNs). The relative num- 
ber densities of the populations and their evolution are con- 
strained by the population fractions in mid-infrared selected 
samples taken from ISO and IRAS (e.g. La Franca et al. 
2004, Matute et al. 2006; see also Manners et al. 2004), and 
are also constrained by early results from Herschel. In mak- 
ing multi-wavelength predictions for each type of object, the 



model uses a set of SED templates, with AGN1 being ap- 
proximately flat in vL u , while the other AGN classes have 
relatively stronger far-infrared excesses. Figure [12] shows the 
detailed predictions of this model for the different classes of 
galaxy, showing how the contribution from the more active 
galaxies is predicted to increase at higher luminosity. As 
with the prediction by Serjeant & Harrison (2005), we again 
find a slight excess in our data at the low-luminosity end 
(although within our la error bars). The dip in the AKARI 
data around L* is not predicted by the model, but occurs 
for both AGNs and SBGs, suggesting it is caused by large- 
scale structure, as discussed above. The shallower slope in 
the AGN population predicted by the model is confirmed, 
supporting the SED assumptions of the model. Neverthe- 
less, our AGN classification is based at optical wavelengths 
so we may in principle be missing AGN2s, for example, and 
could have misclassified AGNs and SFGs if their position on 
the BPT diagram does not reflect their relative bolometric 
contributions. 

In a future paper we will extend our results by present- 
ing additional spectroscopic redshifts in the ADF-S, cover- 
ing more of the field for the mid-infrared-selected popula- 
tion, and also by using Herschel and other data in this field. 
The advent of co-ordinated wide-field far-infrared to sub- 
millimetre surveys, such as AKARI plus BLAST, or Spitzer 
plus Herschel, is enormously extending our ability to con- 
strain the evolving bolometric luminosity function of galax- 
ies. With the additional constraints on AGN dust tori from 
the WISE mission covering both BLAST and Herschel sur- 
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Table 3. The new cluster found in our rcdshift area of ADF-S, 
centred on (04 42 12.5, -53 30 46), within a radius of 45 arcminutcs 
of this centre, and within a redshift range 0.056 < z < 0.061. A 
near position search on the NASA Extragalactic Database (NED) 
yields 649 objects within this area of the sky, although virtually 
all (646) without known rcdshifts. For comparison, other searches 
on NED for the same sized area centred on positions offset by 2 
degrees in RA and/or with reversed Dec typically yielded about 
150-200 sources. 
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vey fields it will be possible to see whether the differences 
we have found compared to prediction are repeated, and to 
improve our assessment of the relative contributions of star 
formation and AGNs to the evolving bolometric infrared lu- 
minosity function. 
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Table 4. Tabulated luminosity function. 
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